Role of bulk and surface phonons in the decay of metal surface states 
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We present a comprehensive theoretical investigation of the electron-phonon contribution to the 
lifetime broadening of the surface states on Cu(lll) and Ag(lll), in comparison with high-resolution 
photoemission results. The calculations, including electron and phonon states of the bulk and the 
surface, resolve the relative importance of the Rayleigh mode, being dominant for the lifetime 
at small hole binding energies. Including the electron-electron interaction, the theoretical results 
are in excellent agreement with the measured binding energy and temperature dependent lifetime 
broadening. 
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Understanding the temporal evolution of quasi par- 
ticles (electron and holes) on metal surfaces is of 
paramount importance to describe many important phe- 
nomena such as the dynamics of charge and energy trans- 
fer, quantum interference, localization and many others. 
This temporal evolution is characterized by a finite life- 
time, r, which refers to the time the quasi particle retains 
its identity. While the lifetime of an excited electron or 
hole is determined by many-body interactions, namely 
electron-electron (e-e) and electron-phonon (e-p) scatter- 
ing processes, the peak width in an experiment might also 
be influenced by electron-defect scattering on crystal or 
surface imperfections However, it was demonstrated 
in recent STM Q and photoemission experiments that 
these defect contributions can be minimized, making it 
possible to analyze the pure lifetime broadening due to 
the formation of a hole in the sp surface state band in 
the i-gap of the (lll)-surface of noble metals. 

These Shockley-type surface states form a two- 
dimensional (2D) electron gas and the e-e contribution 
to the hole lifetime has been rationalized in terms of a 
dominant contribution from intraband transitions within 
the 2D surface state band, screened by the underlying 
3D bulk electron system, and in terms of interband tran- 
sitions (bulk states surface state) On the other 
hand an appropriate calculation of the e-p contribution 
to the lifetime broadening of surface states is still lacking. 
The present work is an attempt in this direction. 

The strength of the e-p coupling is described by the 
electron mass enhancement parameter A, which is, in 
general, energy and momentum dependent. Many prop- 
erties of metals , such as resistivity, specific heat and 
superconductivity, reflect the e-p coupling and can be ex- 
pressed in terms of the Fermi surface-averaged A-value. 
It also reflects the high temperature behavior of the 
broadening T^p = 2TrXkBT, and the e-p contribution to 
the renormalization of the mass m* = m(l + A). The 



anisotropy of A is well known Q and is revealed in e.g. 
cyclotron resonance measurements Q . 

Typically, the phonon contribution to the decay of sur- 
face states is estimated using the Debye phonon model. 
Within this model the Eliashberg spectral function of 
the e-p interaction is proportional to the quadratic den- 
sity of phonon states a'^F{oj) = X{uj/u!d)'^, where ujjj is 
the Debye energy, A is usually obtained from measure- 
ments or theoretical calculations of bulk properties 
However, it is not obvious that this approach should be 
adequate for surface state electrons or holes when the sur- 
face state itself, as well as the surface phonon modes, are 
not taken into account and furthermore, the low tem- 
perature Fep will depend on the model used. A more 
rigorous treatment of the e-p contribution is needed es- 
pecially for surface states close to the Fermi level, be- 
cause for these states the e-e contribution is small and 
the electron-phonon interaction becomes dominant even 
at low temperatures. 

In this letter we present an analysis of the e-p coupling 
contribution to the lifetime broadening of surface elec- 
tron states taking into account all electron and phonon 
states involved in the electron-phonon scattering pro- 
cess, in comparison to new energy and temperature de- 
pendent high-resolution photoemission data. The the- 
oretical analysis is based on a calculation of the full 
Eliashberg spectral function. With this approach we are 
able to resolve in details the contributions from differ- 
ent phonon modes, in particular from the Rayleigh sur- 
face mode and bulk phonons and the general tempera- 
ture dependence. We also obtain the high temperature 
behavior represented by A, which is given by the first re- 
ciprocal moment of the Eliashberg function. We show 
that an approach based on (i) Thomas-Fermi screened 
Ashcroft ion-electron potentials, (ii) one-electron states 
with first principles quality and (iii) a simple model 
phonon calculation give results in good agreement with 
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recently published experimental data ||, ^, Q and data 
presented here. 

The phonon induced lifetime broadening of a surface 
band state with momentum ki and energy cui is given by 

rep(wo) = 2^ / { a^Fj:, (w) [1 + 2 n{Lu) + 
Jo 

f{ujo+uj) - f{uJo - uj)] }duj , (1) 
where is the maximum phonon frequency and 

a'Fj:^ico) = ^ \gf;^\'S{co ~ c.,>) - e^,^), (2) 

where the sum is over final electron states k / of band / 
and phonon modes v with momentum q. 

In Eq. ( H ) we also make use of the quasi elastic 
scattering approximation. The e-p coupling constant is 
given by 

4:;-/^^x(/iE%^(4)-v,/,i.), (3) 

V ^ 

summing over the layers /j of the slab. M and N, are the 
atomic mass and number of ions in each layer of the slab, 
respectively, eq'^^^Rfj.) denotes the complex phonon polar- 
ization vectors normalized over the slab and is the 
screened ion-electron potential. The electron states, of 
the form ip^ ^(a;, z) = {z\n) exp(— fc • a;) , where k and x is 
the electron momentum and electron coordinate parallel 
to the surface, respectively are derived from a slab model 
calculation. This model potential scheme by Chulkov et 
al. [|[ H produces wave functions in good agreement 
with those obtained from ab initio calculations and en- 
ergy spectra reproducing experimental data. 

The Umklapp processes can be neglected as the Fermi 
momentum of the noble metal surface states are small (< 
0.12 a.u.) in comparison with half the minimum recip- 
rocal vector (< 0.75 a.u.). Furthermore, the maximum 
Fermi momentum vector of the final bulk states is < 0.74 
a.u., which implies that the final states are all confined 
to the first Brillouin zone. Thus, the contribution from 
Umklapp processes is less than 10% for small binding en- 
ergies and vanishes completely approaching the F-point. 
To obtain the phonon dispersion relations and polariza- 
tion vectors we perform a calculation using the same force 
constant model as I of Black et al. but in addition, 
we expand the dynamic matrix in Gottleib polynomials 
I p^ in order to optimize the representation of the surface 
phonon modes. The force constant is fitted to reproduce 
elastic constants and bulk phonon frequencies, with sur- 
prisingly good agreement with He scattering experiments 

The screened ion-electron potential is determined by 
the static dielectric function and the bare pseudo poten- 
tial, Vg^{z) = Jdz' i-'^{z,z';q) V^a^S,z'\q) , where q is 



the modulus of the phonon momentum wave vector par- 
allel to the surface and VJ,are 2D-Fourier transform 
parallel to the surface of the bare ion-electron Ashcroft 
pseudo potentials |l^ . We have investigated the quality 
of the screening using the dielectric function according 
to Thomas-Fermi and within RPA (constructed from the 
eigen wave functions and eigen energies of a 31-slab cal- 
culation). The difference is about 1% for both A and Fep 
due to a compensating effect. Referring to the surface 
layer, Thomas Fermi screening is symmetric while RPA 
yields a screening slightly stronger below and slightly 
weaker above. 

We want to point out that in almost all investigations 
of the e-p interaction, the relevant electron scattering 
takes place close to the Fermi surface and only the mean 
e-p coupling, averaged over the Fermi surface, comes into 
play. When considering the lifetimes of surface states, the 
situation is quite different. The e-p coupling becomes 
state dependent as the probed initial electron state is 
fixed. For the surface states of the studied noble metal 
surfaces Cu(lll) and Ag(lll), not only the Fermi sur- 
face is of importance because the binding energy of the 
initial hole state ranges from zero (Fermi level) to about 
0.4 eV. 

The intraband electron scattering (/=«) in Eq.(^) can 
be neglected with the following argument. In the case 
of intraband transitions in the long wavelength limit, 
the matrix element of the gradient of the screened ion- 
electron potential is approximately the expectation value 
of the force acting on the ions in the direction perpen- 
dicular to the surface (V'j ^(x, -z)Vj^ V^{x, z)^. z)) 

K, K, {i\dV^ {z) / dz \i)z. It is well known that 
within linear response, the sum over the forces acting 
on the ions induced by the electron (hole) must be zero 
|]l5| , [l^ , [TtI , |8| , thus the matrix element of the sum of gra- 
dients of the screened ion-electron must be equal to zero 
_F)j=0. Furthermore, the ion displacements forming 
the phonon modes in the surface region, associated with 
small q are locally rigid with a coherence length of ^ir/q 
isotropically in all directions, in parallel and perpendic- 
ular to the surface ||l9| From Eqs. (|^) — (^), it can be 
seen that this restriction and the reduced phase space in 
the intraband scattering process makes this contribution 
negligible compared to interband contribution. 

From Eq. || we see that the Eliashberg function is given 
by the phonon density of states weighted by the electron- 
phonon coupling strength g. In Fig 1. the calculated 
phonon dispersion and the Eliashberg function is pre- 
sented for Cu(lll). The Rayleigh surface mode is split 
of from the bulk phonon band, which gives a lower en- 
ergy peak in the Eliashberg function, at about ~13 meV 
in Cu(lll). The oscillations in the Eliashberg function 
reflects the finite number of layers of the model potential 
(31 layers) in the calculation of electron wave functions 
and thus have no physical significance. 
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The high-resolution photoemission experiments pre- 
sented here were performed with an energy and angular 
resolution of AE w 3 meV and A9 « 0.3°, respectively, 
at VUV photon energy of hv = 21.23 eV (He I). A de- 
tailed description of experimental conditions and sample 
preparation can be found in Ref. ^. A complete data 
set at one temperature, consisting of the photoemission 
data over an angular range of ±6° off normal, was mea- 
sured within the first 15 min after the last annealing step 
at ~ 600° C. As an example, the inset in Fig. 3 demon- 
strates the temperature dependence of the normal emis- 
sion (r point) spectra on Cu(lll) at three different tem- 
peratures. For the quantitative analysis we fitted the 
individual spectra of a temperature series as described in 
Ref. m The resulting Lorentzian lifetime widths F (full 
widths at half maximum: FWHM) are given as open 
symbols in Fig. 3. The binding energy dependence of the 
linewidth at low temperature (cf. Fig. 2) was determined 
from MDC (momentum distribution curve) cuts at the 
respective energy and are given relative to the linewidth 
at the Fermi level Ep, which is nearly completely deter- 
mined by the finite angular resolution. 

The most interesting hole binding energy region is of 
course very close to the Fermi level, in particular when 
the binding energy is less than the maximum phonon 
frequency ujm- In this region we expect that the life- 
time broadening is completely determined by the the e-p 
coupling. The contribution from e-e interaction is very 
small, for Cu(lll) and Ag(lll), Tee < 0.2 meV. 

In Fig. 2 we present the calculated Tep, at T = 30K, 
for Cu(lll) and Ag(lll) together with the experimental 
results. From a simple Debye model we would expect 
a cubic binding energy dependence in the region below 
LUrm which obviously is not the case. The saturation of 
Tep for binding energies exceeding uim is clearly seen in 
experimental data. Adding the contribution from the e-e 
interaction, values close to the experiment are obtained 
(see Tab. |). We note from Fig. 2, that the contribu- 
tion from only the Rayleigh mode gives about 38% of 
Tep beyond the maximum phonon frequencies, indicat- 
ing that bulk phonons give most of the contributions in 
this range. But for binding energies below the maximum 
of the Rayleigh mode energy, this mode alone represent 
on average about the 85% of Tep- 

The main signature of the e-p contribution to the life- 
time broadening is the temperature dependence. For 
binding energies exceeding ksT we can neglect the tem- 
perature dependence of the e-e scattering. The tempera- 
ture dependence of Tep was calculated for the hole state 
in the f point for Cu(lll) and Ag(lll). For a compar- 
ison with the experimental temperature dependence we 
show in Fig. 3 the sum of the constant contribution from 
Tee and the T-dependent Tep. The calculated broaden- 
ings are in good agreement with the experimental data 
points, except for the highest temperature data point for 
Cu(lll). 



In the Tab. ^ we summarize the calculated and mea- 
sured low temperature results for the lifetime broadening 
for the surface states of Cu(lll) and Ag(lll). Summing 
up the calculated contribution from e-e interaction Tee 
and the results of the present work for the e-p in- 
teraction Ten show the proper trend in comparison with 
both STM [| and ARUPS |] data. 

In summary we have presented results from a theo- 
retical and experimental study of the inherent phonon 
induced lifetime broadening of surface states of the no- 
ble metal surfaces Cu(lll) and Ag(lll). Taking into 
account the electron and phonon states, including sur- 
face states is crucial for the understanding of experimen- 
tal findings in particular for states close to the Fermi 
level. The Rayleigh surface phonons are shown to give 
an important contribution to the phonon-induced life- 
time broadening of surface states, in particular for small 
binding energies. We have demonstrated that for noble 
metal surface states, simple Thomas Fermi screening of 
Ashcroft pseudo potentials gives very similar results as a 
more rigorous treatment of the screening. We find excel- 
lent agreements with experiments both with respect to 
the binding energy and the temperature dependence of 
the lifetime broadening. 
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FIG. 3: Lifetime broadening of the Cu(lll) and Ag(lll) sur- 
face hole states as a function of temperature (sohd hue), Pee 
(dotted hue) photoemission data (open circles). The inset 
shows the energy distribution curves (EDC) on Cu(lll) for 
selected temperatures. 

Tables 



FIG. 1: (a) The phonon dispersion from a 31 layer slab calcu- 
lation in the f-M direction of the S.B.Z. (b) The Ehashberg 
function of the hole state in the P point (solid line) and the 
contribution from the Rayleigh mode to the Eliashberg func- 
tion (dashed line). 

FIG. 2: (a) Lifetime broadening of the Cu(lll) surface hole 
state as a function of binding energy. Pee + Pep (solid line). 
Pep (dotted line), the Rayleigh mode contribution to Pep 
(dashed line) and photoemission data (diamonds), (b) The 
same as in (a) for Ag(lll). 



TABLE L Calculated and measured low-temperature lifetime 
widths P at band minimum in meV and mass enhancement 
factors A for the surface states of Cu(lll) and Ag(lll). Pee 
is the calculated lifetime width due to electron-electron scat- 
tering Pep from electron-phonon scattering (present cal- 
culation), FsTM from STM measurements at T=4.6 K 
^ARUPS from ARUPS measurements 

Pep Pee" Pee + Pep PsTA/ ^ ARUPS A 

Cu(lll) 6.6 14.0 20.0 24 23 ± 1.0 0.16 

Ag(lll) 3.6 2.0 5.5 6 6 ± 0.5 0.12 

"i'his value tor i'ee, is a lower limit, because the used model docs 
not take into account the confinement of d electrons |20l . 
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